Anisotropic thermodynamics of d-wave superconductors in the vortex state. 



On 

o 

Q 



o 
o 

di- 
ss 

o3 



i 

G 

O 

o 



> 

O 

m 

On 
O 
OO 
On 
-i— > 



C 

o 
o 



X 



I. Vekhter 1 , P. J. Hirschfeld 2 , J. P. Carbotte 3 , and E. J. Nicol 1 

1 Department of Physics, University of Guelph, Guelph, Ontario NIG 2W1, Canada 
2 Department of Physics, University of Florida, Gainesville, FL 32611, USA 
i Department of Physics and Astronomy, McMaster University, Hamilton, Ontario L8S 4M1, Canada 

(February 7, 2008) 

We show that the density of states and the thermodynamic properties of a 2D d-wave supercon- 
ductor in the vortex state with applied magnetic field H in the plane depend on the angle between 
H and the order parameter nodes. Within a semiclassical treatment of the extended quasiparticle 
states, we obtain fourfold oscillations of the specific heat, measurement of which provides a simple 
probe of gap symmetry. The frequency dependence of the density of states and the temperature de- 
pendence of thermodynamic properties obey different power laws for field in the nodal and anti-nodal 
direction. The fourfold pattern is changed to twofold when orthorhombicity is considered. 



The experimental data accumulated over the last few 
years have established a consensus that the supercon- 
ducting state of the hole-doped high-T c cuprates has a 
predominantly rf-wavc order parameter [Q. Such an or- 
der parameter possesses lines of nodes, which results in 
a gapless excitation spectrum along certain directions in 
momentum space. An important consequence is that the 
low temperature dependence of thermal and transport 
properties of the superconducting cuprates is given by 
power laws, rather than exponential functions with an 
activation energy 

The properties of the vortex state of d-wave super- 
conductors also differ significantly from those of s-wave 
materials: while for s-wave the density of states (DOS) 
and the entropy are determined at low magnetic fields 
H <C H c2 by the localized states in the vortex cores, in 
superconductors with lines of nodes they are dominated 
by the extended quasiparticle states, which exist in the 
bulk along the nodal directions in momentum space pi. 
On the basis of this observation Kiibert and Hirschfeld y] 
suggested a method of calculating thermal and transport 
coefficients in the vortex state microscopically by consid- 
ering only the contribution of the extended states and 
accounting for the effect of the magnetic field on these 
states semiclassically, via a Doppler shift of the quasi- 
particle energy due to the circulating supercurrents. For 
the field applied perpendicular to the superconducting 
layers, the supercurrents can be approximated by the cir- 
cular velocity field around a single vortex, v s = [3/2mr, 
where r is the distance from the center of the vortex, 
and we have set h = 1. Here [3 is a unit vector along the 
current and we use (3 as the vortex winding angle. This 
expression is valid outside the vortex core and up to a 
cutoff of order min{i?, A}, where 2R = 2a _1 y/ <&q/itH is 
the intervortex spacing, $o is the flux quantum, a is a 
geometric constant, and A is the penetration depth. Un- 
der these assumptions the energy of a quasiparticle with 
momentum k is shifted by <5wk(r) = v s • k, and the cal- 
culated physical quantities depend on position and have 
to be averaged over a unit cell of the vortex lattice. The 
results obtained within this framework 0T7|1 describe re- 



cent experiments well [|]-|l3), suggesting that, while the 
effects left outside of this approach are important for a 
quantitative analysis, the method proposed in Ref. Q| 
is qualitatively correct and can be used to analyze the 
properties of the vortex state of unconventional super- 
conductors. 

In this paper we generalize the approach of Ref. Q] 
to consider the experimental arrangement with the mag- 
netic field H in the Cu02 plane, and calculate the density 
of states. We find that it exhibits fourfold oscillations as 
a function of the direction of the applied field, and that 
its energy dependence is drastically different for the field 
directed along the node and anti-node. 

Since the c-axis coherence length, £ c , is shorter than 
the interlayer distance, s, the structure of the vortex state 
for H in the plane differs from that with H|jc, and is com- 
monly modeled by treating the incoherent c-axis trans- 
port as Josephson tunneling between the layers. |l4| |l^] 
The conclusion of Ref. || that the extended states dom- 
inate the thermodynamic properties of the vortex state 
with lines of nodes in the gap remains valid for any orien- 
tation of the magnetic field. The superfluid velocity, v s , 
away from the core is governed solely by the 27r-winding 
of the phase of the order parameter around each vor- 
tex, and, at distances large compared to the size of the 
core, is virtually identical to that of an Abrikosov vortex. 
|I^ , ^5| Since the core size is larger than £ c , the velocity 
field can be approximated by the superflow around a sin- 
gle vortex only for fields H <C Hq = &o/js 2 , where 7 
is the anisotropy ratio; above that field the cores begin 
to overlap. ff5|| While for extremely anisotropic cuprates, 
such as B^S^CaC^Os family, the crossover field is of 
order of a few Tesla, for less anisotropic materials, such 
as YBa2Cu307_<5, Ho > 50 Tesla, and for all experimen- 
tally relevant fields the individual vortices can be treated 
as Abrikosov vortices in the calculation of bulk proper- 
ties. Furthermore, in the regime H <C Hq the intervor- 
tex distance and the structure of the vortex state are 
asymptotically close to those of an Abrikosov vortex lat- 
tice |l6| , suggesting the the approach of Ref. [Q can be 
directly applied to the geometry with the field in the 
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plane. Finally, the c-axis transport remains incoherent 
at low temperatures |l7), and therefore only the ener- 
gies of the quasiparticles with momenta in the plane are 
relevant to the thermodynamic properties and should be 
Dopplcr shifted. 

We now follow the approach of Ref . j| in neglecting the 
contribution of the core states and assuming a spatially 
uniform order parameter Ak over a cylindrical Fermi sur- 
face. In most of this work we consider a pure d-wave an- 
gular dependence of the gap, Ak = A o /(0) = A cos 2(f). 
We consider a magnetic field H applied in the a-b plane, 
at an angle a to the x-axis, and account for its effect on 
the extended quasiparticle states by the Doppler energy 
shift (5ojk(r) = v s • k. The superfluid velocity is approx- 
imated by the flow field of an isolated vortex, which is 
elliptical due to the anisotropy of the penetration depth. 
We can however rescale the c-axis to make both v s and 
the intervortex spacing 2R isotropic, in the London the- 
ory this rescaling is z' = z(X a b/X c ), and since the Fermi 
surface is two-dimensional there is no associated rescal- 
ing of momentum. Then, approximating the unit cell of 
the vortex lattice by a circle of radius R, we obtain 



<5u>k(r) = — — sin/3sin(0 — a), 
P 



(1) 



where p = r/R and E H is the energy scale associated 
with the Doppler shift 



E H = 



(2) 



I-kH 
2~ V^T' 

Here a is a geometric constant of order unity, and in 
the London theory the rescaled Fermi velocity v* = 
VfWXab/Xc, where Vf is the Fermi velocity in the plane. 
In a more general approach v* should be treated as a 
parameter related both to Vf and the anisotropy of the 
vortex lattice. 



The main difference between geometric arrangements 
with the field applied in the plane and that applied along 
the c-axis is clearly seen from Eq. ([[]). For the field applied 
perpendicular to the layers the momentum and real space 
degrees of freedom decouple Q , and the average Dopplcr 
shift is the same at all points on the Fermi surface. In 
contrast, for the field in the plane the average Doppler 
shift becomes dependent on the position at the Fermi sur- 
face, and is given by Eh sin(0 — a). Since quasiparticles 
contribute to the density of states when their Dopplcr 
shifted energy exceeds the local energy gap, an immedi- 
ate conclusion is that the density of states depends sen- 
sitively on the angle between the applied field and the 
direction of nodes of the order parameter. 

To analyze this dependence quantitatively we employ 
the single-particle Green's function which is obtained by 
introducing the Doppler shift into a BCS Green's func- 
tion 101 



G(k, uj n ;r) 



(iuj n - v s k)r + A k Tl + CkT 3 



zv 5 .k) 2 + CI + AI 



(3) 



where u„ is the fermionic Matsubara frequency, £k is the 
band energy measured with respect to the Fermi level, 
and Tj are Pauli matrices. Standard many-body tech- 
niques can be used to compute physical quantities F(r) 
at a fixed position r in real space, and the measured 
quantities are obtained by averaging over the unit cell of 
the vortex lattice according to 



(F)h = - [ pdp f 

7T JO Jo 



df3F(p,P). 



(4) 



We first consider the density of states at the Fermi level 
which is easily accessible experimentally via specific heat 
measurements and is given by 



N (a) = --Um^(TrG(k, w = 0;r)) H 
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sin/3sin(</> — a) 



^/sin 2 /3sin 2 (^ - a) - (A /E H ) 2 p 2 f 2 ((j>). 



(5) 
(6) 



The integrals over p and f3 can be done analytically, and the numerical evaluation of Eq.(||) is trivial. It has been 
shown Q that the nodal approximation, which takes advantage of the fact that the density of states is dominated by 
the contribution of the regions of the Fermi surface near the gap nodes to replace v s ■ k by the Doppler shift at the 
nodes, v s • k„, provides a remarkably good agreement with the numerical results for T, Eh -C Aq. Here it yields 



N {a) 



Eh \ ^ 

A 7T ^ 

nodes 



|sin(0„ - a)\ 



2y/2E H 

A 7T 



max(| sinaj, | cosaj). 



(7) 
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for the field along an anti-node. 
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FIG. 1. The density of states at the Fermi level, No(a), for 
Eh / = 0.1. The solid curve is the full numerical evaluation 
of (5) and the dashed, almost indistinguishable, is the nodal 
approximation Eqn.(6). 

This result was first obtained by Volovik [Q. In Fig.l we 
show the results of full numerical evaluation of the resid- 
ual density of states from Eq.(||) along with the results 
of Eq. (Q) . The density of states exhibits fourfold oscilla- 
tions as a function of the angle of the applied field. There 
is a broad maximum for the field applied in the anti- 
nodal direction, and a sharp minimum for the field along 
the node. Both can be understood if we notice that the 
Doppler shift is given by E\ — Eh \ sin(-7r/4 — a)\ at two of 
the near-nodal regions, and by E2 = Eh\ cos(7r/4 — a)\ at 
the other two nodes. When a field is applied in the antin- 
odal direction, E\ = E2 and all four nodes contribute 
equally to the density of states, as shown in Fig. 2a. On 
the other hand, when the field is applied along a nodal di- 
rection, quasiparticles at that node, which travel parallel 
to the field, do not contribute to the DOS; the Doppler 
shift vanishes at these points. Moreover, since Eh <C A , 
the gap grows faster as a function of the angle (j) near the 
node than the Doppler shift over most of the unit cell 
of the vortex lattice, and therefore the quasiparticle con- 
tribution to the density of states is suppressed over the 
whole near-nodal region, see Fig. 2b. For a field not ex- 
actly in the nodal direction there is always a finite region 
of the momentum space where the Doppler shift exceeds 
the local gap, resulting in a contribution to the DOS and 
sharp minima of N (a). 

Consequently, for a field in the nodal direction two of 
the nodes do not contribute to DOS, while the contribu- 
tion of the other two is a factor of \/2 larger than it is 




FIG. 2. The regions contributing to the density of states for 
the anti- nodal (a), and nodal (b) orientation of the magnetic 
field. 

The density of states is therefore reduced by 30%, in 
agreement with the numerical results. To check how ro- 
bust the oscillations are we numerically computed No(a) 
in a 3-dimensional superconductor and found that the 
results remains qualitatively unchanged, although the 
amplitude of the oscillations is reduced to about 8%, 
in agreement with an analytic estimate. This reduction 
results from incomplete suppression of the contribution 
to the DOS from the nodal lines aligned with the field 
since for quasiparticles outside the equatorial plane the 
Doppler shift does not vanish. We suggest that in real- 
istic materials the amplitude of the oscillations is some- 
where between the two estimates and within the experi- 
mental resolution of the specific heat measurements; the 
published results suggest that this amplitude is of order 
0.3 mJ/mol K 2 at H = 0.005F c2| | @. So far, one re- 
ported measurement performed for two orientations of 
the applied field in the plane did not find the predicted 
oscillations ||. However, an estimate shows that for the 
samples used in Ref. ||] the energy scale Eh for H ~ 8T 
is close to the impurity bandwidth, 7, which may have 
resulted in significant smearing of the fourfold pattern. 
We also note that in an orthorhombic system the induced 
s-wave component of the gap would shift the position of 
the DOS minimum away from the 7r/4 direction, and in 
a heavily twinned crystal, such as used in Ref. ||], this 
would result in rapid filling of the minima and significant 
suppression of the amplitude of oscillations. 

While for a clean sample No(a) oc \/~H independently 
of the angular orientation of H, the energy dependence 
of the density of states depends crucially on the direction 
of the field. For w, E H < A , 

N(u J ,a)~(N 1 {tu,a) + N 2 (uj,a))/2, (8) 

where j| 



Ni(u,a) 



1 



(1 + 2.t 2 ) arcsinx + 3a; VT 



if x = uj/Ei > 1; 
if x < 1 



(9) 
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FIG. 3. The energy dependent density of states for differ- 
ent directions of the applied field and Eh/ = 0.1. 

for i = 1.2, and E\ and E2 were denned above. For the 
field in the antinodal direction 



ttA V ' ZE 2 , 



while for the field along a node 



TV (w, tt/4) 



2E 



H 



TrAr 



2A ' 



(10) 



(11) 



see Fig. 3. The frequency dependence of the density of 
states follows different power laws for the field along a 
node or an anti-node, and in the former case the linear 
slope of N(u>, 7r/4) does not depend on the magnetic field. 
Note that the value of N(0, 0) and N(0, tt/4) differ by a 
factor of \[2 as expected and that the slope of the linear 
term in Eq. (10) is just half the value of the zero field 
case as only two nodes contribute. 

Since the frequency dependence of the density of states 
determines the temperature dependence of thermal and 
transport coefficients, our results have profound effects 
on the properties of clean d-wave superconductors. In 
particular, in addition to the fourfold oscillations of the 
linear-T term in the specific heat as a function of the 
direction of H, the temperature dependence of G jT\Jtl 
is linear for the field in a nodal direction, and quadratic 
for H away from the node. The nuclear spin-lattice re- 
laxation time T\T and superfluid density will also ex- 
hibit fourfold oscillations and a linear or quadratic In- 
dependence depending on the direction of the field. How- 
ever, if the latter quantity is inferred from the penetra- 
tion depth measurements, nonlocal effects due to a di- 
verging coherence length Jl9| in the nodal directions may 
be important. 

Any orthorhombic distortion in the system lifts the 
fourfold degeneracy of the maxima. One possible way to 
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FIG. 4. The effect of orthorhombicity included in the gap 
(d + 0.15s) on N (a) for E H /A = 0.1. The solid curve is the 
full numerical evaluation and the dashed line is the result of 
the first equality in Eqn. (6). 

account for such a distortion is to consider a d+qs super- 
conductor, where g <gc 1. Then, even though the position 
of the nodes is shifted insignificantly, the fourfold pat- 
tern is replaced by two pairs of peaks with the amplitude 
ratio I sin(a — 0.5 arccos(— q))\/\ sin(a + 0.5 arccos(— q))\, 
which differs significantly from unity even for relatively 
small q, as shown in Fig. 4. The anisotropy in the den- 
sity of states and the thermodynamic properties however 
remains robust and is a particularly simple experimental 
probe of the nodal structure of the order parameter. 
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